Astronomy &: Astrophysics manuscript no. MS3921 


February 2, 2008 


(DOI: will be inserted by hand later) 





Evolutionary synthesis of galaxies at high spectral 
resolution with the code PEGASE-HR 

Metallicity and age tracers 

D. Le Borgne^'^, B. Rocca-Volmerange^''', P. Prugniel^, A. Langon^, M. Fioc^'^, C. Soubiran^ 
^ Institut d'Astrophysique de Paris, 98 bts, Boulevard Arago, F-75014 Paris, France 

^ Dept. of Astronomy & Astrophysics, University of Toronto, 60 St. George Street, Toronto, ON M5S 3H8, Canada 
^ CRAL-Observatoire de Lyon, 9 av. C. Andre, F-69561 Saint-Genis Laval, France 
■* Observatoire de Strasbourg, 11 rue de I'Universite, 67000 Strasbourg, France 

Observatoire de Bordeaux, BP 89, 33270 Floirac, France 
^ Universite Pierre et Marie Curie, 4 place Jussieu, 75005 Paris, France 
^ Universite Paris-Sud, 91405 Orsay, France 

Received 2 May 2003 / Accepted 23 June 2004 

Abstract. We present PEGASE-HR, a new stellar population synthesis program generating high resolution spectra 
{R = 10 OOP) over the optical range AA = 400- 680 nm. It links the spectro-photometric model of galaxy evolution 
PEGASE.2 iFioc fc Rocca-Volmerang3ll997fl to an updated version of the ELODIE library of stellar spectra 
observed with the 193 cm telescope at the Observatoire de Haute-Provence (Prugniol & Soubiran 2001a). The 
ELODIE star set gives a fairly complete coverage of the Hertzprung-Russell (HR) diagram and makes it possible to 
synthesize populations in the range [Fe/H] — —2 to +0.4. This code is an exceptional tool for exploring signatures 
of metallicity, age, and kinematics. We focus on a detailed study of the sensitivity to age and metallicity of the 
high-resolution stellar absorption lines and of the classical metallic indices proposed until now to solve the age- 
metallicity degeneracy. Validity tests on several stellar lines are performed by comparing our predi ctions for Lick 
indic es to the models of other groups. The comparison with the lower resolution library BaSeL iLeieune et alJ 
Il997l) confirms the quality of the ELODIE library when used for single stellar populations (SSPs) from 10^ to 
2 X 10^" yr. Predictions for the evolved populations of globular clusters and elliptical galaxies are given and 
compared to observational data. Two new high-resolution indices are proposed around the H-y line. They should 
prove useful in the analysis of spectra from the new generation of telescopes and spectrographs. 
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1. Introduction 

The rapidly increasing number of high-quality spectro- 
scopic galaxy surveys makes it necessary to improve the 
galaxy evolution models used for their interpretation. To 
derive the history of the star formation and the chemi- 
cal evolution of a stellar population from its line-of-sight 
integrated spectrum one can either fit spectrophotomet- 
ric indices and colors with a model, or fit directly the 
observed spectrum with a synthetic spectral energy dis- 
tribution (SED). 

Spectrophotometric indices characterize the strengths 
of spectral features that are sensitive to the age or to 
the metallicity of a stellar populatio n, and generally to 
both. The indices ofE ose et alTl)l994|) and the Lick indices 

Send offprint requests to: Damien Le Borgne, e-mail: 
leborgneSiap . f r 



(IWorthevlll99l iKuntschner fc Davie^llQQSt F&ager et alJ 
:200,q!) are the most widely used. Various models allow us to 
predict the evolution of these indices. They are based on 
assumptions for the stellar evolution (evolutionary tracks) 
and for the history of stellar formation (initial mass func- 
tion, star formation rate. . . ). They usually require prelim- 
inary measurements of t he indices in a library of stars. The 
most recent models ('e.g. lThomas et al .1200.'^ take into ac- 
count the non-solar abundances resulting from the differ- 
ent process es of metal e nrichment by using the response 
functions of lTripicco fc B ell ( 1995). Other families of mod- 
els compute the SED, thus allowing the user to measure in- 
dices a yosteriori if the spectral re s olution is sufficie n t (e.g . 
Fioc fc Rocca-Volmerangel IToQi iLeitherer et 
Eisenstein et^ l2003t iBruziial fc Charlotl l200.ll) . The 
PEGASE code is one of the latter, and its latest version 
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PEGASE.2 l)Fioc fc Rocca-VolmeraM3ll999|) is the basis 
for the present study. 

One of the major issues in studyin g stellar popula - 
tions is the age-metallicity degeneracy llWorthevill99l . 
Most indices are sensitive to both metallicity and age: a 
younger age may be confused with a lower metallicity. The 
comparison of indices featuring Balmer lines with indices 
of metallic lines may in principle break the degeneracy. 
But in practice the tests are difficult and give ambiguous 
results, partly because of the contamina tion by nebular - 
emission, in particular for Ha and H/9 llGo nzalej |l993l 
IWorthev fc Ottavianilll997HKuntschner et al.ll2'oOl|i . 

The degeneracy can also be partly lifted either by us- 
ing broad-ba nd colors and indices in an extended rang e of 
wavelengths (|Gorgas et alJll993t IWorthev et"al]ll994|) or, 
more effici ently, by defining indice s at a higher spectral 
resolution l)jones k Worthevlll995t). PEGASE.2, which is 
based on the librarv of lLeieune et all l)l997lll99q) . is well 
suited to the first approach. Until recently, such SED- 
predicting algorithms were unable to reach high spectral 
resolutions because they lacked adequate stellar spectral 
libraries. 

In a pioneering work, IVazdekis k, Arimotd l)l999|) 

demonstrated the potentiality of using SEDs a t higher res 



olution. They used the library of l.TonesI 1 19981) . thus reach- 
ing R = 2000-3000 (the Lick resolution is about 500), and 
were able to define narrower age-sensitive indices around 
H7 (with band- widths of 10-20 A). 

At high resolution, the analysis of a galaxy spectrum 
must take into account another phenomenon: the inter- 
nal kinematics broadens the lines and reduces the ap- 
parent resolution of the spectrum. Traditionally, measure- 
ments of spectrophotomet ric indices are cor rected for ve- 
locity dispersion (see e.g. iGolev et al.lll999^) . but for gi- 
ant elliptical galaxies with typical velocity dispersions cr 
close to 300 kms~^, the correction would introduce un- 
acceptable errors for narrow indices. For dwarf galaxies 
(cr < 60 kms^^) and globular clusters {a < 10 kms^^), 
corrections would remain acceptable for indices defined on 
2 A-wide passbands. 

Synthetic high resolution SEDs are potentially very 
useful for probing the internal kinematics. They may in- 
deed replace the traditional stellar templates and provide 
both a measurement of the line-of-sight velocity distribu- 
tion (LOSV D) and a constraint on the s tellar population. 
Recently, the lVazdekis fc Arimotol lll999i1 mode l has b een 
coupled to the Ca 11 librarv of ICenarro et alJ l(200lh to 
determine t he LOSVD pr ofiles of bulges of spiral galax- 
ies (Falcon-Barroso et al.L2003) . The LOSVDs are deter- 
mined by deconvolving the observed spectrum with a tem- 
plate spectrum resulting from the population synthesis 
program. Using synthesized spectra to study the inter- 
nal kinematics is not straightforward. Indeed, it is first 
necessary to convolve the template with the spectral in- 
strumental response, which is generally neither Gaussian 
nor constant over the whole wavelength range. To do this, 
the template spectrum must have a significantly higher 
resolution than the observed one. 



To address the questions above, we have coupled the 
last version of PEGASE to a library of high resolution 
{R = 10 000) stellar spectra. In this paper we present the 
resulting code, PEGASE-HR, and focus on its predictions 
of line indices that are sensitive to age and metallicity. 

In Sect. 2, we recall the main characteristics of 
PEGASE.2. We also briefiy describe the mai n features of 
the up dated stellar library ELODIE ( Prugnie l fc SoubiranI 
l2001al) as well as the "interpolator", the method used 
to convert the set of ELODIE stellar spectra to a grid 
of spectra with the regularly spaced stellar parameters 
of the BaSeL library. The synthetic high-resolution spec- 
tra of SSPs are then compared to the corresponding low- 
resolution spectra of PEGASE.2. The variations of the 
high resolution spectra around commonly exploited lines 
are illustrated, as a function of both age and metallicity. 
In Sect. 3, the integrated fluxes, colors and line indices 
predicted with PEGASE-HR are validated by compari- 
son with previous works. In Sect. 4, as a further test of 
the model predictions, we compare them to observations 
of globular clusters and elliptical galaxies. A systematic 
search for new narrow indices is described in Sect. 5. 
Perspectives for the new generation of instruments are 
suggested in the conclusion. 



2. PEGASE-HR: the coupling of PEGASE.2 and 
ELODIE 

2.1. PEGASE.2, a spectral evolution code 

The PEGASE.2 code^ is aimed at modeling the spectral 
evolution of galaxies. It is based on stellar evolutionary 
tracks from the "Padova" group, extended to the ther- 
mally pulsating asymptotic giant branch (AGB) and post- 
AGB phases; these tracks cover all the masses, metallic- 
ities and phases of interest for galaxy spectral synthesis. 
For a given evolutionary scenario (typically characterized 
by a star formation law, an initial mass function and, pos- 
sibly, infall or galactic winds), the code consistently com- 
putes the star formation rate and the metallicity of gas 
and stars at any time. The nebular component (contin- 
uum and lines) due to Hii regions is roughly calculated 
and added to the stellar component. Depending on the 
type of galaxy (disk or spheroidal) , the attenuation of the 
spectrum by dust is then computed using the output of a 
radiative transfer code; this code takes into account scat- 
tering (^Fioc 1997). 

PEGASE.2 uses the BaSeL llLeieune et al.ll997lll99?j) 
library of stellar spectra and can therefore synthesize low- 
resolution {R ~ 200) ultraviolet to near-infrared spectra 
of Hubble sequence galaxies as well as of starbursts. In 
PEGASE-HR, the BaSeL library is replaced by a grid of 
spectra interpolated from the high-resolution ELODIE li- 
brary of stellar spectra. 



^ Version 2 of PEGASE, Projet d'Etude des GAlaxies 
par Synthese Evolutive in French. The code is available at 
http: //www. iap . f r/pegase/. 



Le Borgne et al.: Galaxy evolution at high spectral resolution 



3 



2.2. ELODIE, a high resolution library of stellar spectra 

The ELODIE library is a stellar database of 1959 spec- 
tra for 1503 stars, observed with the echelle spectrograph 
ELODIE on the 193 cm telescope at the Observatoire de 
Haute Provence. Previous versions of the library were pre- 
scnted in Soubiran et al. (1998) and Prugnicl & Soubiran 
It has been updated for the present work by 
doubling the number of spectra, which greatly improved 
the coverage of the parameter space (in effective temper- 
ature, surface gravity, and metallicity) . The data reduc- 
tion has also been improved, in particular the flux cali- 
bration, and the wavelength range has been extended to 
AA = 400-680 nm. For the purpose of population synthe- 
sis, the original resolution {R = 42 000) has been reduced 
to i? = 10 000 at A = 550 nm, or more precisely to a 
gaussian instrumental profile of FWHM~ 0.55 A over the 
whole range of wavelengths. The typical signal-to-noise 
(S/N) ratio of the spectra is 500/ A. The connection of sev- 
eral echelle spectra, each one defined in a small wavelength 
range, is needed to build a single spectrum over the entire 
range AA = 400-680 nm. The complex procedure used to 
perform this task leads to an overall (broad-band) photo- 
metric precision of 2%. The precision reaches 0.5% when 
the ELODIE spectra are normalize d to a low resolution 
continu um, as described in detail in lPrugniel fc SoubiranI 
120014. 

The atmospheric parameters of the stars (Toff, g, 
[Fe/H]) ar e taken from the literat u re in an up-to-date 
version of ICavrel de Strobel et al.l l|200l|) . when avail- 
able. Quality weights are applied to the parameters 
when several estimations w- ere available (for details, see 
IPrugniel fc SoubiranI '200 laV. Otherwise they are esti- 
mated with the TGMET procedure, which consists in 
least square fits of the target spectrum to reference spec- 
tra with known atmospheric parameter s llKatz et alJl998t 
ISoubiran et"aLl l2003l) . Unlike Worth ev et all I^1994^ we 
chose not to calibrate the tempera ture of the giant stars 
on their spectral type as given bv iRidgwav et alJ l)l980l) 
because this study is contested in several recent pa p ers 
l|Sekiguchi fc Fukugitallioool: iHoudashelt et al.ll2000(l . As 
a consequence, our estimated T^s for the giant stars is 
globallv smaller by a bout 5% than the values given by 
IWorthev et all l(l994l) . This wiU be discussed in Sect.O 

The HR diagram coverage with the estimated param- 
eters is extensive: 0.27 < logjo(5 [cm.s~^]) < 4.97 and 
3185 K<roff< 55 200 K, with -3.21 <[Fe/H]< 1.62. This 
coverage is illustrated in Fig. 

The up-to-date version of the library is available on- 
line''. Fully reduced spectra, as well as the estimated stel- 
lar parameters, are provided. 



Since the number of sp ectra included in the ELODIE 
library has doubled since IPrugniel fc SoubiranI l)2001al) . 
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^ and also | http: / /www.obs.u-bordeauxl .fr /public/astro/ 1 
CSO /elodie Jibrary.html 
^ http:/ /www-obs. univ-lyonl.fr/~prugniel/soubiran/ 
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Fig. 1. Distribution in the T^s- logj^Q g diagram of the 1503 
stars in the ELODIE library in four bins of metallicity. 



we update here the comparison of the Lick indices mea- 
sured on__tlie_stellar spectra with the reference indices 
from IWorthev et"al. (1994) for the 187 stars in common. 
Figure 121 illustrates this comparison for 16 indices, show- 
ing the generally good agreement between the two sets of 
spectra. However, a closer look reveals a general trend for 
our strongest indices to be weaker than the reference in- 
dices, which is expressed by the slopes being smaller than 
unity. To investigate thi s bias, we also compared our spec- 
tra to the Jones hbrary i Leitherer et aljn"99fiH,ToneJl998l) 
for the 314 stars in common. This library has the advan- 
tage of being well calibrat ed in wavelength, as opposed 
to the Lick/IDS spectra of lWorthev et alJ l|l994l> . When 
compared at the same resolution (1.8 A), the two sets of 
spectra agree remarkably well. 

Table ^ gives the quantitative comparison of these 
three libraries in terms of Lick indices. The indices of the 
Jones library were measured by ourselves when comput- 
ing this table. The good agreement between ELODIE and 
Jones spectra is shown by the slopes being closer to unity 
and the smaller dispersions a. It is worth noticing that the 
indices published by Jones differ from the ones we mea- 
sure from their spectra: they seem to have been artificially 
corrected in order to match the Worthe v et al.. (,1994^1 in- 
dices by changing the overall slope. Our routine used to 
measure Lick indices was tested successfully on seven stars 
provided by G. Worthey. To measure the indices on the 
higher resolution spectra of the ELODIE or the Jones li- 
braries, we use the wavelength-dependent resolution given 
in Wort hev fc Ottaviani (1997) . 

The small bias observed between the measurements 
made on our library and the Jones library is probably due 
to the imperfect calibration of the spectra in either library. 
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Fig. 2. Comparison of Lick indices measured on the ELODIE spectra with the reference Lick/IDS indices 
(Wort hev et al.lll994|) for 187 stars in common. The dotted Unes are y = a; ± ctids where ctids is the uncertainty 
on the Lick/IDS indices. For each panel, we give the slope of the linear regression and the dispersion a around it. 
Ideally, cr/triDS - 1- 



The reason for the discrepancy between our measurements 
and the Lick/IDS ones is still unclear, but it might be 
attributed to a systematic effect in the wavelength scale 
of the Lick library. 

In the following, we do not correct any of our indices 
for this very small bias. Section will show that its con- 
sequences on SSPs are usually small, although they may 
become important when red giants dominate the spectra. 

To use the set of ELODIE spectra as an input stel- 
lar library to PEGASE, the spect ra were interpolated 
on the BaSeL i Leieune et al.lll99^ grid of parameters 



(Tcff, logj^Q [Fe/H]) with pi ecewise polynom ials similar 
to the Lick fitting functions ijWorthevI Il994j) . in various 
regions of the HR diagram. The interpolated spectra were 
then normalized to match the mean BaSeL fluxes in the 
band 5500-6000 A. The detailed description of the inter- 
polator, as well as the comparison of the BaSeL library 
with the interpolated stellar library, are described on a 
web site*^ and will be publ ished in a forthcoming paper 
llPrugniel fc Soubiranll200^ . 



* http:/ /www.obs.u-bordeauxLfr/public/astro/CSO/ 
elodie_library.html 
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Table 1. Comparison of Lick indices measured on 
ELODIE stellar spectra with indices from other databases 
of stellar spectra, fo r stars in common. The indices for 
IWorthev et al.1 were taken from the table published 

in their paper. The indices for the Jones library were mea- 
sured by ourselves (see text). The values are the mean 
offset (ELODIE— other database), the slope of the linear 
regression and the dispersion a around it (see Fig. 
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The range of stellar parameters in the BaScL library 
being much larger than the range of parameters covered by 
the ELODIE library, we chose not to use the extrapolated 
spectra of the interpolated library. To exclude the extrap- 
olated spectra, we set a minimum level for the density of 
real stellar spectra at each point of the BaSeL grid. The 
resulting interpolated library then contains fewer spec- 
tra than the BaSeL library (2690 instead of 4422), but 
each of its spectra is reliable. A minor modification in the 
PEGASE algorithm was necessary to use this new grid. 

2.3. SSPs at high spectral resolution 

Figure 13 presents a prediction of a 10 Gyr-old high reso- 
lution SSP spectrum at solar metallicity. 

The labeled grey areas in the plot iden tify the pass- 
bands of some Lick indices l|Worthev et al ][l994). We can 
notice that each passband includes many lines of various 
chemical elements. More precise line identifications and 
their evolution are presented in the next section. To high- 
light the spectral resolution improvement, we have over- 
plotted the same SSP computed with PEGASE. 2 using 
the BaSeL library. 

It is worth noting that the adopted evolutionary tracks 
are parametrized by their total metallicity Z , while the 
stellar library is labeled by the iron abundance [Fe/H]. In 
the present work, we assume the relation log]^Q(Z/0.02) = 
[Fe/H]. The possible correlation between metallicity and 
relative abundances of a-elements in the spectral library. 




5300 



5380 



Fig. 4. Comparison of a portion of a smoothed PEGASE- 
HR SSP spectrum (Z(7i , 3 Gyr) with the model of 
iBruzual fc Charlod ll2003i) . The IMF is common to both 
models l|Sa]net.eT<ll955t) .' a) Both models overplotted, at 
the same low resolution R ~ 2 000. The flux is given in 
units of 10~^ L0.A~^.Mq-^. b) Relative difference of the 
fluxes at low resolution (i? ~ 2 000) c) The original high 
resolution PEGASE-HR SSP spectrum [R = 10 000). 

due to the selection of stars from the solar neighborhood, 
will be invoked in the following to explain some of the 
differences with other models. 



2.4. The spectra around specific stellar lines 

Figure El presents the normalized spectral distributions of 
SSPs centered on the main stellar lines (H(5, H7, H/9, Mgb, 
FeA5270, FeA5335, NaD). The qualitative effects of age 
and metallicity can be assessed: in the left column, all 
SSPs have solar metallicity and their ages are at 1, 4, 
13 Gyr; in the right column, all SSPs are 10 Gyr-old and 
their metallicities are [Fe/H] = —1.7, —0.4, 0.4. Lick index 
passbands are shown with vertical lines and the grey ar- 
eas correspond to the pseudo-continuum passbands. The 
widths of the Lick passbands contrast with the narrow 
profiles of the stellar lines, outlining the possibility of re- 
fined indices. 

2.5. The library of SSPs 

Further examples of synthetic SSPs spectra at the resolu- 
tions of PEGASE.2 and PEGASE-HR are shown in Fig.lHl 
For clarity, PEGASE.2 spectra are offset vertically. The 
qualitative agreement in the energy distribution and the 
broad features is very satisfactory. A quantitative com- 
parison is provided below. Note that PEGASE-HR and 
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6400 6500 6600 ^ 6700 6800 

Fig. 3. PEGASE-HR high-resolution spectrum of a 10 Gyr-old SSP of solar metallicity compared to the low resolution 
spectrum from PEGASE.2. The flux is given in units of 10"^ Lq.A-^Mq^ where Lq = 3.826 x 10^^ erg.s"^ Grey 
areas correspond to the passbands of some Lick indices. 



PEGASE.2 are based on the same algorithm, so that the 

evolution of global parameters (metallicity, mass of gas, 
stellar mass and all outputs independent of the library) is 
unchanged. 
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Fig. 6. Evolution of main stellar lines and continua for an SSP as a function of age (ages = 1,4, 13 Gyr, from light gray 
to black) at Zq (left) and metallicity ([Fe/H] = —1.7, —0.4, 0.4, from gray to dark) at 10 Gyr (right). Line identifications 
are from the revised ILLSS catalogue (.Coluzzi 1999). Grey areas show the blue and red pseudo-continua of the Lick 
indices and solid vertical lines delimit their central passband. In some cases {H'fA, ^^a), significant evolution with age 
is observed within the pseudo-continua (see text). 
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Fig. 5. Comparison of PEGASE.2 and PEGASE-HR 
spectra of SSPs, at metallicities 2.5xZq and Zq/5, and 
ages 0.5 Gyr and 10 Gyr. The PEGASE.2 spectra are 
plotted with a vertical offset (a constant in each panel) 
for clarity. 



3. Consistency of the PEGASE-HR outputs 

3.1. Quantitative comparison to PEGASE.2: 
luminosities and colors 

Figure [7| presents a comparison of the integrated fluxes 
predicted by PEGASE-HR and by PEGASE.2 through 
the quantity 

/ ^6700A \ 
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for SSPs of various ages (1 Myr to 20 Gyr) and metallic- 
ities (from [Fe/H] = -1.7 to [Fe/H] = 0.7). Absolute dif- 
ferences for ages greater than 10 Myr and [Fe/H]> —1.7 
are always smaller than 0.04 magnitudes. 

Figure |S1 compares the slopes of the energy distribu- 
tions produced by the two versions of the code. The slope 
is measured with a color index, based on filters with rect- 
angular passbands spanning 4300-4600 A ("blue" band) 
and 6400-6700 A ("red" band). Color differences between 
the two models, over the whole range of metallicities and 
ages greater than 10 Myr, are lower than 0.12 magnitudes. 
Once more, the largest disagreements are confined to the 
youngest ages and lowest metallicities. 

Overall, the colors and fluxes of PEGASE-HR are con- 
sistent with those of PEGASE.2 for ages greater than 
10 Myr and [Fe/H]> -1.7. This result confirms that the 
BaSeL and ELODIE libraries are very similar at low res- 
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age (Myr 

Fig. 7. Comparison of the integrated optical flux pre- 
dicted by PEGASE.2 and PEGASE-HR for SSPs using 
SF (4200-6700 A). The size of the symbols scales with 
magnitude difference (see text in Sect. l3.1l for details). 

olution, at least over the range of parameters in which 
stellar spectra contribute significantly to the optical light 
of a galaxy. Such an agreement was not obvious a priori 
since the two libraries are constructed with very differ- 
ent methods: the former is based on theoretical spectra, 
color-corrected to fit observed stellar continua, whereas 
the latter is purely empirical. As a consequence, the low- 
resolution extended SEDs predicted by PEGASE.2 and 
the high-resolution spectra of PEGASE-HR can be used 
together to refine SED studies. 



3.2. Comparison with Lick indices for SSPs 

Before studying lines at high resolution, we first want to 
check that our predictions at medium resolution arc con- 
sistent with those of other models. 

The traditional approa ch for the computat i on of 
Lick indices of SS Ps (e.g IWorthev fc Ottavianil Il997t 
iBressan et aLlllQQ^ is based on a library of indices of in- 
dividual stars and on analytic representations of this set 
of data (the fitting functions). Thanks to the resolution 
of the ELODIE spectra, our approach is more direct: we 
compute the synthetic spectrum of a galaxy or SSP, de- 
grade the spectral resolution to the one of the Lick index 
definitions (8 to 11 A depending on wavelength, as de- 
scribed bv lWorthev fc Ottaviani 1997.) . and then measure 
the indices directly on the smoothed spectrum. This allows 
us to easily predict indices for any complex evolutionary 
scenario. In the following, we compare our results for SSPs 
with those of other groups. 

3.2.1. Metal-sensitive indices 

The left hand side of Fig. |^ shows the variations with 
age of the metal-sensitive Lick indices Fe5270, Fe5335, 
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Fig. 9 . PEGA SE-HR Lick indices (solid lines ) for SSPs compared to models (dotted lines) of IW orthev fc Ottavianil 



(W95V lBressan et al. (199 6) (BCT96). lT"homas et all l)2n03l) (TMB03) and Bru ziial & Ch ariot (2003) (BC03). 
The figures on the right hand side of each panel represent the quantity [Fe/H]. For BC03 only, the dashed line shows 
the indices computed with the Geneva tracks instead of the Padova 1994 tracks, for solar metallicity. 
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Fig. 8. Comparison of the optical colors predicted by 
PEGASE-HR and PEGASE.2 for SSPs. The adopted color 
compares the fluxes integrated over 4300-4600 A in the 
blue and 6400-6700 A in the red. The size of the triangles 
scales with the color difference expressed in magnitudes 
(5color) between the two models. Head-up triangles are 
plotted when PEGASE-HR produces a redder slope than 
PEGASE.2. 



Mgh and NaD as com puted with PE GASE-HR, by 



Bressan et all (Il996l). by iThomas et all |2003 ). and by 

it 



BruzualfcCharlotl ll2003() . The cornpariso ns are shown at 
various metallicities. 'Tho mas et ahl l|2003|) provide results 
for various relative abundances of the g-e lements, a nd we 
have plotted those for [a/Fe] = 0. The ISalneterl 1^19^^ 
initial mass function (IMF) is used in the three models 
(power law index x — —1.35). The lower and upper stellar 
masses are adapted fo r each compari son. 

PEGASE-HR. 'Bressa n et all (|l996!) and 
iBruzual fc Charloil |2003) all use the Padoya st ellar 
tracks in this comparison, while Irhomas et al.l l)2003(l use 
what we will call the "Cassisi tracks" , a compilation from 



Cassisi et all lll997l):lBono et al.l l|l997|) : lMaraston (199^ . 
alasnich et al.1 ll200(ll). Moreov er, both iThomas et all 



2003() and 'Bres san et iaJI 1 1996() use "fitting functions" 
to derive their synthetic indices, contrarily to us and to 
iBruzual fc CharlotI ll2003|) . 

The agreement of the iron indices between the three 
models is quite good, except for the youngest SSPs at the 
lowest metallicities. This shows tha t ELODIE spectra and 
the iron index fitting functions of IWorthev et all l|l994l) 
are compatible over the range of parameters relevant for 
the synthesis of stellar population s, as expected from 
the dire ct comparison performed bv'Prugni el fc SoubiranI 
I 2001a|^ . As stated in Maraston et al... (,2QflJ), the Cassisi 
tracks should lead to slightly lower values of the indices. 
This is what we observe, especially for iron indices. The 
discrepancy with Bruzual fc Chariot (2003) at super-solar 
metallicity will be explained below. 

The situation for the Mgh index is more complex, 
since magnesium is an a-element. Both the ELODIE li- 



brary and the Lick fitting functions rely on an empir- 
ical stellar data set, in which one should expect the 
anti-correlation between [Fe/ H] and fa/Fe] characteriz- 
ing th e solar neigh borhood l|Fiihrmannl 1 19981 iFulbrigliil 
2002). IThomas et al. (2003|) have included a correction 
for the local anti-correlation between [Fe/H] and [a/Fe] in 
their sub-solar computations. At solar metallicity, where 
no such correction is required, their predictions and those 
of PEGASE-HR agree very well. As expected, PEGASE- 
HR produces a higher Mgj, index at sub-solar metallicities. 
We have obser ved that our predictio ns match better the 
values given bv IThomas et all l)2003|) for [a/Fe] =0.3. 

The differences betwee n the Mgf, predictions of 
PEGASE-HR and those of iBressan et all ^996), which 
are significant at super-solar metallicities, can only result 
from selection effects in the underlying stellar samples or 
in peculiarities of the interpolations in regions of the HR 
diagram where only few stars are available. Although the 
differences in stellar samples are probably the main source 
of disagreement between the models at sub-solar metallic- 
ity, the latter option is likely to explain the strong differ- 
ence observed at [Fe/H] — 0.4, as demonstrated below. 

We have compared the Mgf, indices measured on the 
ELODIE spectra (after interpolation onto the grid of fun- 
damental parameters used in the population synthesis 
code) with those obtained for the same grid of parame - 
ters from the fitting functions of IWorthev et al T lll994l) . 
The average offset is negligi ble, as can be expected from 
IPrugniel fc SoubiranI l)2001a|) and Fig. [2 However, the in- 
dices measured on ELODIE spectra of non-giant stars 
(Mgf, < 6) have a slight tendency to have lower values 
than those derived from fitting functions for strong lines, 
and higher values for weak lines. This effect, also present 
for the individual stars (see Fig.jSJ, is of the order of a few 
tenths of an Angstom, i.e. not sufficient to explain the dif- 
fere nces observed f or Mg;, in Fig.Elbctwcen PEGASE-HR 
and lBressan et all ((1996) . More important is that the Mgf, 
fitting functions rise dramatically with decreasing temper- 
atures for low gravity red giant stars (Tcff< 3800 K). When 
using the Padova stellar evolution tracks, these stars con- 
tribute significantly to the optical emission of old metal- 
rich populations . Very few such stars were present in the 
library used by IWorthev et aD ^9^), all of them with 
sub-solar metallicity, and their Mgh indices are high in- 
deed (values above 15 A, up to 18.5 Aj. The Mb;, indices 
of the coolest super-solar giants in the ELODIE library are 
lower by up to several A: the strongest Mg^ index for a gi- 
ant is Mgf, = 14.3 A. Also, there are only two giant stars 
in the ELODIE hbrary with [Fe/H]>0 and Mgh > 7 A 
(HD18191 and HD169931). The interpolation of the stellar 
spectra is very hazardous with so few stars in this range 
of metallicity, effective temperature and surface gravity. 
Therefore, one should not be surprised to observe strong 
differences in the indices of a-elements for old, super-solar 
stellar populations. 

Moreover, a comparison of our estimate d effective tem- 
peratures with the values published in IWorthev et alJ 
((19941 for the stars in common shows that the calibration 



Le Borgne et al.: Galaxy evolution at higli spectral resolution 



11 



in 
in 



3.0 






PEGASE-HR 




2 








BC03 






2.5 






BCT 




1 ^ 












2.0 




.V\V 






- 






'■.\x 




X 


-1 - 


1 .5 




















-2 : 










1 .0 










— 3 





5 


10 15 











age 


(Gyr) 






4.0 




, , , , 






5.5 








_ ^ ' 




5.0 ; 


3.5 










4.5 ; 










X) 


A C\ - 
■M- . U ■ 


3.0 




//,.•' 




Ol 


3.5 - 












2.5 










3.0 - 












2.5 - 


2.0 










9 



5 10 15 
age (Gyr) 



5 10 15 
age (Gyr) 



5 10 15 
age (Gyr) 



Fig. 10. Lick indices for [Fe/H] = 0.4 SSPs The i ndices of 
this work (solid linel. iBruzual fc Charlol l)2003|) (dotted 
line) and iBressan et alJ l)l996f) (dashed line) are shown. 
The shaded region shows the uncertainty in the indices 
associated to the uncertainty in the effective temperature 
of the giant stars (see text for details). 



of Tcs for the cool giant stars is uncertain by about 150 K 
(5%) which could also explain some of the differences be- 
tween the models at super-solar metallicities. The effect of 
this uncertainty on the predicted indices at [Fe/H] = 0.4 
is comparable to the discrepancy between the models: 
Fig. Uni shows four indices of a [Fe/H] = 0.4 SSP ob- 
tained by shifting artificially the temperatures of the stel- 
lar tracks by ±5 % (shaded region). The o t her models at 
this metallicity ijBruzual fc Chariot 120031: iBressan et al.l 
predict values compatible with this uncertainty. 
However, this shaded region probably overestimates the 
error due to the uncertainty in stellar effective tempera- 
tures: the error in T^ft for the dwarf main sequence stars 
is smaller than 5%. At solar or sub-solar metallicities, the 
contribution of the giant stars to the spectrum of an SSP 
is smaller. The uncertainty in Toff then becomes negligible 
in comparison to the uncertainty in the stellar tracks. 

The index NaD is in quite good agreement with the 
other models, which is remarkable because this index 
is contaminated by many telluric lines. We took great 
care to remove these lines from the high resolution stel- 
lar spectra. The non-monotonic evolution of the NaD in- 
dex with metallicity for old SSPs in the predictions of 
iBruzual fc CharlotI ((200.ll) reflects the difficulty to model 
this index correctly. 

3.2.2. Hydrogenic indices, age tracers 

Balmer lines are generally considered to be excellent age 
tracers. Ha and H/3 are often strong absorption lines, but 



they may be filled in by circumstellar and interstellar emis- 
sion. Thus, they provide less reliable age diagnostics than 

H7 and H(5 which are much less aff ected. 

To analyze stellar populations, IWorthev fc Ottavianil 
defined Lick indices around easily measurable 
Balmer lines: H/?a, H'Ja and US a for A stars with broad- 
ened wings, and iifip, H'jp and RSp for later types, 
with narro w-er Balmer lines. T hey u se the Revis e d Yale 
isochrones ijCreen et all Il987l) and IVandenberd 
isochrones extrapolated, when necessary, to some parts 
of the HR diagr am. The IMF adopted by them has 
a ISalpete slope, with Mniin — 0.1 M© and 

Afmax = 2 Mq (they model old populations only). 
Figure IHl (right hand side) shows an overall satisfactory 
agreem e nt between PE GASE-HR. Wor thev fc Ottaviani 
" 19971) . iThomas et all (|2003) and .Bruzual fc Chariot 
20031) for metallicities [Fe/H]> -1.00. However, we are 



also in rather good agreement with IBruzual fc Chariot 
(,2003.) for [Fe/H] = -1.7. We beheve that some of 
the discrepancies in the predictions of Lick indices for 
Balmer lines at low metallicity are due to different treat- 
ments of the helium -bu rning stars in the stellar tracks 
l|Charlot et al.lll996t lYil 12005^ . This is illustrated in the 
right panels of Fig. [HI where the IBruzual fc Charloil l(2003l) 
predictions for solar metallicity are g iven with the Padova 
tracks, but also with the Geneva ijSchaller et al.l 119921: 
ICharbonnel et aLlll996() tracks. As for metal-sensitive in- 
dices, many of the discrepancies in the predictions of the 
Balmer Lick indices, in particular at super-solar metal- 
licities, are likely due to differences in the stars used to 
build the stellar libraries and in their estimated effective 
temperatures. 



4. Synthesis of evolved populations 

4.1. Globular clusters 

Predictions of spectral feature strengths in SSP spectra 
can be directly compared to observational data of co- 
eval stellar populations, as found i n globular clus t ers. Th e 
data used here wer e obt ained bvlBica fc AlloinI 1 1986aj) . 
iHuchra et all l)l996(l and lTraeer et all l)l998(l " 



Equivalent width of the H/3 line 

Equivalent widths (EW) of the H/3 line were measured as 
a function of age and metallicity for a sample of 63 ob- 
ser ved star clusters in the G alaxy, the LMC, and the SMC 
bv iBica fc Alioh] lll986albl) . We predict the evolution of 
the H/3 EW with PEGASE-HR, with the same spectral 
resolution (12 A, Gaussian smoothing) and the same EW 
definition: 



EW 



[l-i^(A)/F,(A)]dA 



w 



(1) 



where Fc is the continuum in the middle of the win- 
dow W. The continuum is the upper envelope of the 
spectrum: maximal fluxes around the points AA — 
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Fig. 11. Evolution sequences with age of the H/3 equiv- 
alent width from PEGASE-HR for various metallicities 

i solid lines) compared t o the binned globular cluster data 
Bica fc Alloh]ll986alb^ . 

4570, 5340, 6630 A (±20 A) are connected by straight 
lines. EW is given in A. 

Figure ^2 plots the binned data compared to the 
PEGASE-HR sequences of H/3 equivalent widths for vari- 
ous metaUicities ([Fe/H] = -1.7. -0.7, -0 4.0., 4). Note 
that the isochrones used bv lBica fc AlloinI (1986b) to esti- 
mate ages differ from the ones used here, allowing no very 
precise age comparison. 

The models agree with the data at ages roughly be- 
tween 1 and 10 Gyr. At younger ages, emission lines might 
partly fill in absorption lines in the observations, and on- 
going star formation may rule out the model of the strictly 
instantaneous (<1 Myr) burst described by an SSP. At the 
oldest ages, the metal deficiency of globular clusters is ex- 
treme, and we reach the regime in which our interpolated 
spectral library suffers from larger uncertainties. 



The iron-H/? diagram 

Among the 13 i ndices measured in 1 93 extragalactic glob- 
ular clusters bv lHuchra et al.l l)l996|) . we select the index 
Fe5270 as a representative metallicity tracer and H/3 as 
the best available ag e tracer (see also Ta ble Et. The def- 
initions of indices in iHuchra et alJ l)l996l) differ from the 
Lick ones: 



-2.5 log 



2Fj 



{Fci + Fc2) 



(2) 



where i^/, Fci and Fc2 are mean fluxes in the pass- 
band, blue and red continua respectively. The passbands 
are also different from the Lick ones. In Fig. ^1 the 
grid of PEGASE-HR models (with ages=0.3-14 Gyr and 
[Fe/H] = —1.7 to 0.4) is compared to the data sample. 
The grid encompasses a wide range of data, when taking 
into account the large observational error bar. However, 
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Fig. 12. Fe 5270 vs. H/3 indices f or 193 extragalactic globu- 
lar clusters l|Huchra et al.ll996() . Sohd hues are PEGASE- 
HR SSP models for several ages and metallicities, as la- 
beled on the sequences. The typical error bar for the data 
is plotted above the grid. For M33 globular clusters, the 
size of the full squa res indicates [Fe/H] (as determined by 
iHuchra et alJ 119961) . following the [Fe/H] scale given by 
the empty squares plotted on the 14 Gyr isochrone. 

several points (mostly M31 clusters) are outside the the- 
or etical grid. Assumin g that the behavior of the indices 
of IHuchra et al.l l)l996l) is comparable to that of the cor- 
responding Lick indices, the offset in H/3 for the low- 
est metallicity clusters is consistent with the uncertain- 
ties already mentioned for this regime. For M33 clusters, 
Brodic & Huchra (1991) give estimations of the individual 
metallicities. Overall, we find a rather good agreement be- 
tween measured metallicities of clusters and the grid. This 
agreement is illustrated by the size of the squares in the 
plot (filled squares for M33 data and empty squares for 
the model predictions). 

More recently, some precise measurements of Lick in- 
dices were made on globular c lusters of the Milky way 
and M31 and were compiled bv lTrager et al.l l)l998() . The 
measurement error is smaller for these data, and enables 
a more precise determination of ages. Figure El shows the 
observed values of Fe5270 and H/3 together with our syn- 
thetic grid for these indices. It appears that the locus of 
the globular clusters in the grid is approximately the same 
as for the extragalactic clusters of Fig. E] However, the 
ages are determined much more precisely and they span a 
narrower range (3-20 Gyr) with mainly sub-solar metal- 
licities. 



4.2. Elliptical galaxies 

Elliptical galaxies are often modeled as purely coeval stel- 
lar populations since the bulk of their stellar population 
formed within a very short timescale (less than 1 Gyr) . We 
hereafter prefer to use the typical scenarios of the model 
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Fig. 13. Fe5270 vs. H/3 ind ices for 50 globular clusters of 
the Milky Way aiid M31 l|Trager et al.lll998(l . The grid 
corresponds to PEGASE-HR SSP models, with iso-age 
sequences (dotted lines) and iso-[Fe/H] sequences (solid 
lines) . The median error bar for the data is plotted below 
the grid. 



PEGASE l|Fioc k Rocca-Volmerangelll997i) including in- 
fall and galactic winds. The star formation rate is pro- 
portional to the mass of gas, with respective timescales of 
1 Myr for starburst, less than 1 Gyr for elliptical and up 
to 10 Gyr for spiral galaxies. 

Figure El shows the predictions of the main 
Lick indices with our three scenarios (spiral, el- 
liptical, and instantaneous b urst, as defined in 
iLe Borgne fc Rocca-Volmerangd l2002j) compared to 
the same indices meas ured on local elliptical galax- 
ies l|Trager et al ] 1200(1 . The assumed galaxy age is 
13±0.5 Gyr, but the indices are not very sensitive to age 
in this range. Mgf, predictions are systematically low, 
which is consistent with the known over-abundance of 
a-elements in elliptical galaxies. 



5. A systematic search for new narrow indices 

Lick indices are considered, to a certain extent, to be good 
estimators for age (H/3 index) or metallicity (Fe or Mgb 
indices). Yet these indices are defined for medium resolu- 
tion (FWHM ~ 8 A): many absorption lines, from various 
elements, are included in each passband. The full resolu- 
tion power of PEGASE-HR makes it possible to build new 
narrow indices for a better estimation of ages and metal- 
licities. At the limit, indices could be as narrow as a few 
tenths of an Angstrom, but they will be accessible only 
to a few instruments. Also, the exposure time would have 
to be very high to obtain a high S/N ratio. Moreover, 
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Fig. 14. Lick indices of PEGASE-HR models for eUiptical, 
Sc and single burst g alaxies together wi th data from lo- 
cal elliptical galaxies jTrager et al.ll2000(l . The age of the 
data is arbitrarily set by us to 13 ± 0.5 Gyr, randomly 
distributed (for the clarity of the plot). The presence of 
emission lines associated to ongoing star formation ex- 
plains the low values of the H/3 index predicted for evolved 
spiral galaxies. 
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stellar velocity dispersion broadens the lines (50 kms^^ 
corresponding to FWHM=2.0 A at 5000 A for example). 
Extremely narrow indices would turn out to be unusable. 
For these reasons, we restrain our search to 2 A- wide in- 
dices. In order to disentangle age and metallicity, we tried 
to build couples of indices sensitive to either age or metal- 
licity, producing the most orthogonal and regular grids 
in the age-Z plane, with the smallest measurement errors 
at a fixed r nean S/N ratio (here after 50). Indeed, as em- 
phasized bv lCardiel et al.l l|200,ll) . this association of small 
measurement errors (which depend directly on the defini- 
tions of the indices) with an orthogonal grid is the most 
promising combination for disentangling age and metal- 
licity signatures. The errors in in dex measurements ar e 
estimated following the formulae of lCardiel et al.l l)l998l) . 

5.1. Evolution of the pseudo-continuum 

The definition of the continuum is an important aspect 
of a nv index defini t ion. T he continuum can be defined, 
as in lBica fc AlloirJ l|l986a|^ . by the upper envelope of the 
spectrum, but is then strongly sensitive to noise and emis- 
sion lines. The Lick indices define a "pseudo-continuum" : 
a straight line between the mean fluxes in a blue and a 
red window. This choice is more robust and is convenient 
for the analysis of observations. The Lick normalization 
is empirical because the pseudo-continuum may evolve by 
itself with age and metallicity of the underlying popula- 
tions. Fig.llSlshows this effect on the H7 line. SSP spectra 
at various metallicities are presented at 10 Gyr with the 
Lick index resolution. Two normalizations are presented: 

a) at 4448 A (top panel); the blue part of the Lick pseudo- 
continuum dramatically changes because it falls in the 
G band; 

b) at the Lick pseudo-continuum (bottom panel); the 
variations seem to reside in the H7 line itself. 

Depending on normalization, conclusions may strongly 
differ. 



5.2. A systematic search method 

We have developed a method to find narrow indices, de- 
fined like the Lick indices, that could be good estimators 
of age or of metallicity. The first step is to compute SSP 
spectra for various metallicities, and various ages. Then, 
we normalize each spectrum to the maximum value within 
a sliding 50 A-wide window. We thus obtain spectra with 
almost flat continua. Next, we measure, for each wave- 
length, the relative flux variations induced by either age 
changes or metallicity changes. We use the relative dif- 
ferences between these variations to select narrow feature 
windows and pseudo-continuum windows; we thereby de- 
fine indices particularly sensitive to age or to metallicity. 
Since the S/N ratio of the stellar library is finite, we take 
into account the uncertainties in the measurement of the 
flux variations relative to age or to metallicity. This leads 
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Fig. 15. Evolution of the II7 line as a function of metallic- 
ity for a 10 Gyr-old starburst: a) normalization at 4448A; 
b) normalization to the Lick H'ja pseudo-continuum. In 
both panels, the spectrum is degraded to the Lick reso- 
lution. Grey zones represent the 117^1 Lick blue and red 
continuum definitions and solid vertical lines delimit the 
117^ feature passband. 

us to select easily measurable indices only (i.e. with small 
measurement error bars). 

5.3. The age-sensitive index H^_VHR 

The H7 line, as a Balmer line, is a good age-sensitivity 
tracer. Moreover, unlike Ha and II/3, it is seldom filled 
in by nebular emission (H7/H/3 ~ 0.5 in Hii re- 
gions). The upper panel of Table [21 presents various 
definitions for a H7 index (wavelengths and resolu- 
tions): H7_F (Worthey. Lick resolution ~ 8 A), Il7_Vaz 
l|Vazdekis fc Arimoto h99^. a = 200 kms^i; FWHM ~ 
6.8 A), H7_HR (J ones fc Worthev (1995), ct = 80 kms-^: 
FWHM ~ 2.7 A). We hereafter propose the new index 
H7_VHR (VHR standing for Very High Resolution) with 
FWHM=2.0 A (cr -60 kms"!), weU suited to PEGASE- 
HR. 

Figure [TEl presents the comparison of the age-Z grids 
built up with various H7 indices for the age axis, and 
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Table 2. Various H7 indices and resolutions (last column). H7_VHR and H7_Z are the new ones that we propose (see 
Sect. 10 and [Oil. 



Name 


Feature 


Blue band 


Red band 


FWHM (A) 


H7_F 

H7JHR 

H7_Vaz 


4332.500 4353.500 
4338.607 4342.347 
4332.000 4352.250 


4284.750 4321.000 
4333.000 4335.000 
4331.000 4340.750 


4356.000 4386.000 
4348.000 4350.000 
4359.250 4368.750 


9.5 
2.7 
6.8 


H7_VHR 
H7_Z 


4337.600 4341.200 
4351.000 4354.000 


4327.600 4336.400 
4318.000 4322.500 


4348.800 4366.800 
4362.000 4365.000 


2.0 
2.0 



with the Lick <Fe> index for the [Fe/H] axis. The last 
grid (H7_VHR vs. <Fe>) shows that the H7_VHR index 
is particularly insensitive to the metallicity and allows a 
slightly better age estimation than H7_HR at low metallic- 
ity. For old, solar-metaUicity stellar populations, it is very 
similar to the H7_HR index, but the error bar associated 
to its measurement (at S/N=50 on the plot) is smaller. It 
may be measured at high (> 50) S/N ratio with the new 
generation of instruments. 



5.4. A metal-sensitive index: Hj_Z 

The high spectral resolution (FWHM=2.0 A) makes 
it possible to present the new H7_Z index, highly sen- 
sitive to metallicity . This mi ght look surprising, but 
IVazdekis fc Arimotol l)l999l) and ljones fc WorthevI (l995|) 
already noted that many metallic lines are present in 
every Lick index. We call this new index HjJZ because 
it is very close to the H7 line, but the central passband 
actually measures a Chromium (Cr) line. This index is 
voluntarily defined for a spectral resolution corresponding 
to FWHM=2.0 A. This narrow index is quite comple- 
mentary to the index H7_VHR: the passbands are very 
close to each other, and a single short spectrum is enough 
to perform an age-Z diagnostic. The independence of the 
parameters traced by the two indices is exemplified in 
Fig. El The diagram H7_VHR (age-sensitive) vs. H7_Z 
(Z-sensitive) should be useful to interpret high-resolution 
spectra. 



5.5. Z and age sensitivities 

We present in Table|3|updated values of the sensitivities to 
age and metallicity of some Lick indices (H7, 11/3, Fe5335, 
Fe5270) and of the old and new H7 indices. We use the def- 
inition of Worthey (1994.) for the quantity (Aago/Apc/H]) 
called "Z-sensitivity" . It is the ratio of the fractional 
change in age to the fractional change in Z needed to in- 
duce the same variation of the index value. A value close 
to zero means that the index is completely insensitive to 
metallicity variations, whereas a large absolute value in- 
dicates that the index is insensitive to age variations. We 
also give these fractional changes Aago and Apc/u] for a 
10% variation of the index value around the zero point 
(Zq, 12 Gyr). A comparison to values obtained with pre- 



Table 3. Sensitivity of line equivalent widths to age (1- 
13 Gyr) and metallicity ([Fe/H] = -1.7 to 0.4). A[fo/h] 
represents the fractional change in [Fe/H] needed to pro- 
duce a 10 % variation of the index value at the zero 
point (Z0, 12 Gyr). Aagc is the similar quantity computed 



with the age. 



i[Fe/H] 



is the ratio of fractional change in 



age to the fractional change in [Fe/H] r equired to pro- 
duce the same change i n the index value l)Worthevlll994 
I.Tones fc Worthevlll995^ . High absolute values correspond 
to high sensitivity to Z. This quantity is evaluated around 
the (Zq, 12 Gyr) zero point. The last c olumn reproduces 
the values in Table 6 of I WorthevI l)l994 and in Table 2 of 
iJones fc WorthevI lll99,'C 



Index 


A[Fe/H] 


^agc 


Aago 


f Aago \ 


A[Fo/H] 


\^iF=/Hiy w,j 


Fe5335 


0.41 


0.75 


1.8 


2.8 


Fe5270 


0.42 


0.74 


1.8 


2.3 


H7_Z 


0.73 


1.03 


1.4 




Mgt 


0.55 


0.56 


1.0 


1.7 


BSa 


0.16 


0.14 


0.9 


1.1 


H7A 


0.36 


0.31 


0.9 


1.0 


H7F 


0.23 


0.15 


0.7 




H7_Vaz 


-3.04 


-0.42 


0.1 




H/3 


-1.62 


-0.26 


0.2 


0.6 


H7_HR 


6.50 


-1.06 


-0.2 


0.0 


H7_VHR 


11.05 


-1.06 


-0.1 





vious models l)Worthevlll994t Ijones fc WorthevI Il995|) is 
also given in the last column. Our new index H7_VHR is 
quite insensitive to metallicity variations around the zero 
point. The new index H7_Z is very sensitive to variations 
of metallicity, even more than the Mgj, index. 

6. Discussion and Conclusion 

The code PEGASE-HR is a spectro-photomctric evolu- 
tion model of galaxies with very high spectral resolution 
in the optical. It inherits the PEGASE.2 accurate mod- 
eling of galaxy evolution, as demonstrated by the com- 
patibility of colors and low-resolution spectra between the 
two models. The ELODIE library is complete enough to 
give an extensive coverage of the HR diagram. Efforts are 
currently made to complete the library with low metallic- 
ity stars. Its quality depends on the estimation of stellar 
parameters with TGMET and on the interpolator which 
will be soon published. The ELODIE archive continues to 
grow from new observations, and the density of stars in 
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:[Fe/H] = -1.7, 




Fig. 16. H7 indices vs. <Fe> Lick index. From top to bot- 
tom, the H7 indices are computed with the definitions of 
Worthey, Jones, Vazdekis, and our Very High Resolution 
index (see text for d etails). Error bars are computed from 
ICardiel et all l|l998>) with an average S/N ratio=50. Solid 
lines are iso-age sequences (1, 5, 12, 17 Gyr). Dotted lines 
are iso-[Fe/H] sequences (-1.7, -0.7, -0.4, 0., 0.4). The 
thickness of the lines increases with age and [Fe/H]. Notice 
the particularly flat iso-age sequences of the H7_VHR in- 
dex, which indicates a very low Z-sensitivity. 




Fig. 17. H7_VHR index (age tracer) vs. HjJZ index 
(metallicity tracer) indices. Error bars are computed from 
ICardiel et alJ l)l998|) with an average S/N ratio of 50. Same 
symbols as in Fig. El 



the parameter space {Tcs, log^Q^j Pe/H]) increases. The 
improvement will be particularly noticeable for fast evolv- 
ing stars and a better accuracy will be obtained for young 
starbursts (age < 10 Myr). 

We checked that PEGASE-HR, degraded at low res- 
olution, is in good agreement with other recent models: 
P EGASE.2. Bressan ct al. (1996), Thomas et al. (2Qq3), 
or iBruzual &: CharlotI l|2003j) . The predicted evolution of 
Lick indices for Balmer lines and metallic lines is compat- 
ible with most other models and observations. The small 
number of cool giant stars with super-solar metallicity in 
the ELODIE library and the uncertainty on their effective 
temperature makes the predictions of indices quite uncer- 
tain for old, metal-rich stellar populations. The influence 
of the non-solar enrichment for a-elements in the ELODIE 
library needs more investigation. 

We systematically explored the wavelength domain to 
find new high resolution indices, sensitive to either age or 
metallicity. Inside the passband of the classical II7 Lick 
index we find two very high resolution indices: one is sen- 
sitive to age only , and is quite similar to the index of 
Ijones fc WorthevI l|l995tl . The other one is sensitive to 
metallicity only. These 2 A-wide features are defined for 
velocity dispersions up to 60 kms~^. We note that these 
indices, despite a systematic investigation at high resolu- 
tion, are not fundamentally different from the previously 
existing ones. This may indicate that future improvement 
will mainly come from high-resolution SED fitting rather 
than from classical indices. 

This model is suited to the analysis of high resolution 
spectroscopic observations with rest-frame wavelengths 
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falling in the 4000-6800 A interval. In particular, the 
most recent spectrographs on large telescopes (GIRAFFE, 
VIMOS or ISAAC at the VLT, CMOS or GNIRS at 
the Gemini Observatory, EMACS on Magellan, TWIN at 
Calar Alto, ISIS at WHT. . . ) will benefit from this model. 
We already presented tests of inversion methods show- 
ing their ability to separate disk and bulge components 
from a spectrum, deriving simultaneousl y their popula- 
tion c haracteristics and their kinematics I Prugniel et alJ 
l200.ll lOcvirk et all 120031) . This approach will be funda- 
mental in many scientific applications, such as the study 
of bulges for w hich correcting the d isc contamination is 
a pre-requisite l|Prugniel et al.ll200l|) . or of high surface 
brightness dwarf objects. 

The code, as well as SSPs, can be downloaded on 
the PEGASE web site at http://www.iap.fr/pegase/. 
Because the output spectra now contain ten times as many 
wavelength points as the PEGASE. 2 spectra, the SSPs are 
provided in FITS format. They contain all the information 
previously given by PEGASE. 2 outputs, except the colors. 
Measurements of Lick indices are included. A code to con- 
vert the FITS files into the PEGASE.2 ASCII format is 
also available for the user's convenience. 

We stress that this new tool is unique for studying 
stellar populations in nearby and distant galaxies ob- 
served with the new high resolution spectrographs. Other 
improvements will come from multi-wavelength analy- 
ses, by combining PEGASE-HR optical spectra with far- 
ultraviolet and near-infrared SEDs. Another improvement 
will come from the implementation of non-solar abun- 
dances in the evolution. For this purpose, we are presently 
measuring chemical abundances from the high resolution 
spectra of the ELODIE library. 
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